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1. 


INTRODUCTION 


In  this  report,  a  coax-fed  monopole  antenna  completely  immersed 
in  water  is  investigated  theoretically  and  the  results  are  corroborated 
experimentally.  The  major  thrust  of  the  investigation  is  to  acquire 
a  better  understanding  of  the  properties  of  the  current  on  an  antenna 
in  a  lossy  medium.  Coupled  integral  equations  are  derived  for  the 
current  on  the  antenna  and  the  electric  field  in  the  aperture  formed 
where  the  coaxial  guide  joins  the  conducting  ground  plane  (Fig.  1). 

These  equations  accurately  characterized  the  above-mentioned  current 
and  electric  field  for  the  structure  illustrated  in  the  figure. 

Numerical /analytical  techniques  needed  to  solve  thse  equations  are 
developed  in  the  report.  Numerical  results  are  obtained  and  are 
compared  with  experimental  data  available  elsewhere.  These  equations 
and  the  numerical  techniques  developed  for  solving  them  provide  the 
basis  upon  which  one  may  found  an  accurate  analysis  of  an  antenna  in 
a  lossy  medium. 

The  monopole  antenna  illustrated  in  Fig.  1  is  the  extension  of  the 
center  conductor  of  a  coaxial  waveguide  whose  outer  conductor  terminates 
at  the  conducting  ground  plane.  The  monopole  is  taken  to  be  a  conducting 
tube  (open-ended)  with  a  vanishingly  thin  wall,  it  is  fed  or  driven  by  a 
time-harmonic  signal  (e^1")  from  the  coaxial  guide  and  receives  its 
excitation  through  the  annular  aperture  where  the  coax  opens  into  the 
ground  plane.  The  monopole  is  of  height  h  and  radius  fi  and  it  resides 
in  a  semi-infinite,  homogeneous  medium  characterized  by  (vi,e)  with 
e  (■  e0e  -j  — )  allowed  to  be  complex  to  account  for  losses.  The 


radii  of  the  inner  and  outer  coax  conductors  are  a  and  b,  respectively, 
and  the  medium  in  the  annular  region  between  the  conductors  may  be 


uniform  in  which  case  it  is  characterized  by  (yc,£c)  or  it  may  be  a 

gas  (y  ,£  )  in  which  a  dielectric  (y  ,£  )  bead  of  length  L  is  inserted 
8  8  c  c 

to  provide  mechanical  support.  (Gee  Fig.  2  .)  In  the  former  case  of 
a  uniform  solid  dielectric,  L  is  made  to  approach  infinity.  As 
suggested  in  Fig.  2  ,  the  monopole  radius  may  differ  from  that  of  the 
coax  inner  conductor. 

We  wish  to  calculate  the  current  induced  in  the  monopole  and  to 
determine  the  load,  or  admittance,  which  the  monopole  presents  to  the 
coaxial  line.  The  known  excitation  is  specified  to  be  a  TEM  wave 
traveling  in  the  positive  z  direction  in  the  coax,  whose  field 
components  are 


E 


i 

P 


h: 


a) 


where  k2  *oj2y  e  and  n2  ■  1J  /e  and  where  En  is  the  known  complex  ampli- 

tude  of  the  TEM  electric  field  incident  upon  the  annular  aperture.  If 

the  coax  is  gas-filled  and  the  center  conductor  is  supported  by  a 

dielectric  bead,  k  and  ti  are  replaced  by  k  and  n  ,  respectively.  The 

8  8 

viewpoint  taken  here  is  that  the  monopole  can  be  treated  as  a  scatterer 
which  is  illuminated  by  the  field  radiated  by  the  coaxial  aperture  in 
the  screen.  This  approach  enables  us  to  apply  aperture  and  scattering 
theory  to  the  present  analysis  and  provides  the  basis  for  the  formulation 
of  a  pair  of  coupled  integral  equations  for  the  monopole-coax  structure. 
As  an  aid  in  the  derivation  of  equations,  the  formulation  is  partitioned 
into  three  parts:  field  in  coaxial  region,  field  radiated  by  annular 
aperture,  field  scattered  by  monopole. 


Cross-section  of  coax-fed  monopole  above  a  ground 
plane  (for  coax  filled  with  uniform  solid  dielectric 
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2.  FORMULATION  OF  INTEGRAL  EQUATIONS 

Coupled  integral  equations  are  formulated  in  this  section  for  the 

unknowns  I  and  Ea  where  I  is  the  total  axial  current  on  the  monopole  and 

P 

Ea  is  the  p  component  of  the  electric  field  in  the  annular  aperture  A 
through  which  the  coaxial  guide  and  monopole  antenna  are  coupled.  These 
integral  equations  are  derived  by  requiring  (-t)  that  the  <j>  component 
of  the  H-field  in  the  coax  be  equal  to  that  in  the  exterior  region  when 
both  are  evaluated  in  the  limit  as  the  points  of  observation  appr  ,h 
the  annular  aperture  and  ( -Lt)  that  the  electric  field  tangential  t  r.he 
monopole  and  to  the  conducting  plane  be  zero.  For  this  purpose  t 

cl 

magnetic  field  is  expressed  in  the  coaxial  guide  as  a  function  of 
and  of  the  excitation  and  in  the  exterior  region  as  a  function  of 

Ea  and  I. 

P 


Magnetic  Field  in  Coaxial  Guide 


.co 


The  magnetic  field  in  the  coaxial  guide  can  be  written  in  terms 
of  the  coax  excitation  of  (1)  and  aperture  field  E  as 


H^°(P,z)  -  Hac(p,z)  +  H^(p,z) 


(2a) 


where 


and  in  which 


H*(p,z) 


Ep (p ’ )GC (p , z ; p ' )dp ' 


Hr<p,z)  “  2  n^P  cos  V 

c 


-  • 


(2b) 


(3a) 


-»■-«  •  * 
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is  the  so-called  short-circuit  magnetic  field  which  would  exist  in  the 
uniform-dielectric  coax  if  the  aperture  at  z  =  0  between  the  coaxial 
guide  and  the  exterior  region  were  short-circuited.  If  the  coax  is 
gas-filled  with  a  bead-supported  center  conductor,  the  short-circuit 
magnetic  field  in  the  region  -L < z < 0  would  be 

r,  jk  L 

sc  ^0  e  ® 

H,  (p,z)  =  2  —  - 5 — — — ■ — : - : — r — =-  cos  k  z,  z  e(-L,0).  (3b) 

d>  ’  p  n  cos  k  L  +  i  n  sm  k  L  c 

g  c  J  c  c 

In  the  case  of  the  gas-filled  coax,  (2a)  is  valid  only  in  the  bead 
material  (-L  <  z  <  0) ,  and  it  is  assumed  that  higher  order  modes  excited 
at  the  discontinuity  at  z  =  0  are  vanishing  small  at  z  =  -L.  The 
Green's  function  GC  is  determined  by  standard  methods  [  1  ]  to  be 


ik  z  -j2k  (z+L) 

_c.  ..  1  ^  c  1-Te 

G  (p,z;p  )  *  -  .  e  T 


n  p&n- 

c  a 


1  +  r  e 


-j  2k  L 
c 


where 


where 


—  p’  l  7-p-f-r  *  <P’>  ^  *<P)  e3XqZ 

nc  q-l  XqWq  dp  q  uP  q 

ze(-L,0) 


-  n 


no  +  n  * 
8  c 


V0)  -  VV0  W’  -  Jo(V)  VV’  • 


and  where 


,  k2  <  k2 
tq  c 


-jv'F  -k2  ,  k2  >  k2 

J  tq  c  tq  c 


6 


in  which  k^,  q“l»2,...»  is  the  root  of  $Q(b)»0.  In  (6),  and 

N  are  the  vC^  order  Bessel  and  Neumann  functions  and  the  norm  N2  in 
V  q 

(4)  is 


Nq  *  |  ^  4,to)]  dC‘T- 


_2_ 

TT2  ‘ 


(8) 


The  Green's  function  Gc  of  (4)  is  valid  for  both  the  uniform-dielectric 
coax  and  the  gas-filled  coax  with  bead-supported  center  conductor;  in 
the  former  case,  F  is  set  equal  to  zero  in  the  first  term  of  (4). 

Field  in  Exterior  Region 

The  magnetic  field  in  the  exterior  region  is  determined  as  the  sum 


H®X  =  H*"  +  H® 

<P  <P  <P 


(9) 


where  is  that  radiated  by  the  annular  aperture  in  the  absence  of  the 

g 

monopole  and  is  that  scattered  by  the  monopole.  The  field  incident 
upon  the  monopole  is  that  radiated  by  the  aperture.  As  an  aid  in 
deriving  expressions  for  the  exterior  field,  the  exterior-field 
equivalent  model  of  Fig.  3  is  introduced.  In  this  figure  one  sees 
that  the  annular  aperture  A  at  z  =  0  is  short-circuited  and  that  an 

A 

equivalent  surface  magnetic  current  of  density  M  =  M^<p  where  = 

or  =  -E®  is  placed  on  the  shorted  annulus  [  2 ] .  This 
model  and  its  field  are  equivalent  to  the  original  structure  and  its 
field  in  the  exterior  (z>0)  region  [  2  ].  By  image  theory  one  arrives 
at  the  final  model  of  Fig.  3 ,  which  is  used  as  a  guide  in  exterior 
formulations. 


7 


Will'll. 


(a) 


(b) 


Lg.  3.  Exterior-region  equivalent  models  of  coax-fed  monopole. 
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In  the  absence  of  the  monopole,  the  field  (Er,Hr)  due  to  the 


magnetic  current  is  determined  from 


r  1 
E  =  -  —  V  x  F 

£  -  - 


(10a) 


Hr  -  -j  jjy  (k2F  +  TV  •  F) 


where 


(10b) 


*J1 


.  -Jkl£-P’ I 

2M,  (p')<j>*  - -  dS' 

*  '  r-p'l 


or  where  F  ■  F,$  with 

<P“ 


e  f  a  f  g-J^l  z  ■ 

F  (p,z)  -  -f-  p’E®(p')  cos*1  - - 

9  2  P  L  fz2  + 


-jkfz2  +  R2 ] 2 


dcfr'dp’ 


[z  +  R2] 


(12a) 


in  which  k2  ■  u)2ye  and  n2  *  y/e  and 


R2  *  [p2  +  p,z  -  Zpp’coscj)']  . 


(12b) 


From  (10)  and  (12)  one  readily  finds  that  the  field  components  Er 

and  are 
<P 


_  1  S  f  a  e-JKlz  ■ 

E*<p>2)  ■  1?  p  J  p'EP(p,)  J  cosp’  — Z77 


.2  .  D  2  i  *5 


-jk[z2  +  R2] 


d<f>'dp' 


[z2  +  R2] 


2  _i_  1,21*5 


r  -jk[z  +  R2] 

hI(P,z)  -  j  -T-  -  p'EP(p')  cosp’  - - r -  d$'dp\  (1*) 

<•  2"  n  J  P  J  [z2  +  r!]>! 


The  field  (E  ,H°)  due  to  the  total  axial  current  i  induced  on  the 


(15a) 


ES  =  -j  (k2A  +  W  •  A) 

and 

Hs  «  I  Vx  a  (15b) 

M  “  - 

where  A  has  only  a  z-component  given  by 


VP’Z)  -  47 


TT  -jk[(z-z’)2  +  R2] 


21*5 


-h 


-7T  [(z-Z1)2  +  R2] 


,21*5 


d<f>'dz' 


(16) 


in  which  R2  is  R2  of  (12b)  with  p'  replaced  by  A.  As  is  consistent  with  the 

a. 

fact  that  the  excitation  of  the  monopole  emanating  from  the  aperture 
is  <|>  independent,  the  current  on  the  monopole  and  the  resulting 


scattered  field  are  also  <p  independent.  Expressions  for  the  field 
s  s 

components  E ^  and  follow  immediately: 


E*(p,z)  -  -j  — 


4TTk 


-h 


IT  -jk[(z-z’)2  +  R2] 
e _ 

^  [(z-z’)2  +  Rj]*5 


2-1  *5 


(17) 


d<J>'dz' 


tt  -jk[(z-z')2  +  R2  ]** 


J  x<.->  ^  J 

— h  — TT  'l 


d<{)  'dz ' 


(18) 


Integral  Equations 

The  final  integral  equations  for  E^  and  I  result  from  enforcement 
of  boundary  and  transition  conditions.  The  boundary  condition  is  that 
the  total  electric  field  tangential  to  the  monopole  and  ground  plane 
surfaces  be  zero.  The  tangential  E-field  on  the  ground  plane  is 
zero  by  virtue  of  the  use  of  image  theory  in  the  construction  of 


V.  1 


expressions  for  _E  and  El  and  that  on  the  surface  of  the  tubular 
monopole  is  made  zero  by  enforcement  of  the  following  equation: 

Es(*,z)  +  Er(/L,z)  -  0  ,  ze(-h,h)  .  (19) 

z  z 

The  required  transition  condition  is  that  the  magnetic  field  tangential 
to  the  annular  aperture  A  in  the  exterior  region  and  that  in  the  coax 
be  equal  when  the  points  of  observation  approach  a  common  point  in  A. 
This  condition  is  ensured  upon  enforcement  of 


lim  H*T°(p,z) 
ztO  41 


py 

lim  H.  (p,z)  ,  peA 

z+0  ^ 


or  of 


lim  H*7°(p,z) 
z+0  41 


lim 

z4-0 


H*(p,z) 


H*(p,z) 

; 


ieA 


(20) 


Eqs.  (19)  and  (20)  must  hold  simultaneously  and,  upon  utilization  of 
(2),  (3),  and  (13)— (18) ,  they  lead  to  the  coupled  integral  equations 

below  for  E  and  I: 

P 


-j 


4irk 


fr  +  k2 

dz2 


I(z')  G“(A.,z;z')dz' 


4 


-h 

b  -i 

► 

fe’J 

P'Ep(P')  Ga(p,z;p')dp' 

L  a  J 

■  0  ,  ze(-h,h)  (21a) 
p-4 


and 


_ 1_  9_ 

4tt  3p 


I(z')  G  (p,0  ;z')dz'  +  j  £ 


-h 


n  J 


P’E“(p')  G  (p,0;p,)dp' 


E*(p’)  GC(p,0;p’)dp'  -  H®C(p,0),  pe (a,b) 


(21b) 


11 


in  which 


TT  -jk[  (z-z')2  +  R^]** 


G  (p,z;z') 


-TT  [<«-«'>*  +»J] 


^  d$' 


(22a) 


and 


-jk[z2  +  R2]*5 


Ga(p,z;p’) 


COS<f>' 


-TT 


[  z  2  +  Rz  ] 


TT**' 


<22b) 


For  subsequent  convenience,  the  integral  equations  are  written  below 
in  operator  form: 

E®[I;z]  +  E^[Ea;z]  -  0  ,  ze(-h,h)  (23a) 

z  z  p 

and 

H®[I;p]  +  H*[E*;p]  -  Hj[Ea;p]  -  H*c(p,0)  ,  pe(a,b)  .  (23b) 

I  and  Ea(— M^)can  be  determined  by  solving  the  coupled  integral  equations 
and  from  knowledge  of  these  two  quantities  the  radiated  field  and  the 
input  TEM  reflection  coefficient  of  the  coaxial  guide  can  be  computed. 


3. 


NUMERICAL  SOLUTION  METHOD 


£ 

The  integral  equations  for  Ep  and  I  are  far  too  complex  to  be 

solved  by  any  but  a  well-conceived  numerical  method.  Due  to  the 
c 

presence  of  G  ,  one  must  sum  an  infinite  series  comprising  terms 
involving  derivatives  of  which  depends  upon  Bessel  and  Neumann 
functions  whose  arguments  are  determined  by  a  solution  of  a  transcen¬ 
dental  equation  for  each  q.  In  addition  the  following  terms  in  the 

integral  equations  involve  integrals  whose  integrands  may  be  singular: 

s  r  c  gx 

E  (,K,  z) ,  E  C't.z),  lim  H,(p,z),  and  lim  H,‘  (p,z).  A  technique  has 

2  2  z+0  *  z+O  ♦ 

been  developed  for  solving  numerically  the  integral  equations  and 

schemes  to  handle  all  the  difficulties  enumerated  have  been  devised. 


This  technique  is  described  below. 

The  current  I  on  the  monopole  and  its  image  as  well  as  pE  in 


the  annulus  are  represented  by  linear  combinations  of  pulses  11^  of  the 


forms 


Kz)  -  l  i  n  (z) 
n-1  n  n 


and 


N 

a 


where  I  and  V  are  unknown  constants  to  be  determined,  N  is  the 
n  n  s 

number  of  non-zero  pulses  in  the  interval  (-h,h),and  N  is  the  number 

cl 

of  pulses  in  the  interval  (a,b).  The  pulse  functions  are  defined  as 


1,  5  e  (Cn-Ac/2,  Cn  +  AC/2) 

- 
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0,  otherwise 


where  C  is  the  pulse  center  and  A£  is  its  width.  On  (-h,h)  the 
n 

n*'*1  pulse  center  is  at  z  ■  -h  +  nAz  and  the  pulse  width  Az  has  value 

Az  *  2h/(Ng+ 1)  while  on  (a,b)  the  nC^  pulse  center  is  at  «  a  +  Ap(n-*s) 

and  the  pulse  width  Ap  has  value  Ap  ■  (b-a)/N  .  In  keeping  with  the 

boundary  condition  l(±h)«0,  "half  pulses"  are  placed  on  the  sub  intervals 

(-h,-h+Az/2)  and  (h-Az/2,h),  and  their  coefficients  are  set  equal  to 

zero.  Eq.  (21a)  is  tested  [3]  with  triangles  A  ,  m«l,,..,N  ,  defined 

m  s 

by 


In  (27b)  K(z-z')  =  GS(4.,z-z')  is  the  exact  kernel  [  4  ]  of  cylindrical 
antenna  and  scatterer  theory.  Eq.  (21b)  is  tested  with  delta  functions 

or  is  simply  enforced  at  match  point  Pm  located  at  the  centers  of  pulses 
to  arrive  at 


N 

s 

l  I 

n=l 


j,as 

n  mn 


N 

a 

+  l  V  [Y  -YC  ] 

w  n  nm  h>»>  •* 

n=l 


H 


sc 

m 


(28a) 


where  HSC 
m 


m-  1,2,. ...N 

a 

sc 

H<j)  (pm’0)s  where  p0  *  a  +  Ap(®Js),  and  where 


ras 

mn 


n 

Bp  j  V*’’  GSC0.0;2-)dzn 

-h  JP  = 


m 


(28b) 


mn 


(p’>  Ga(pm,0;p')dp' 


(28c) 


and 


mn 


U 

j  nn(p*)  ±  Gc(pm,0;p')dp’  . 


(28d) 


In  view  of  the  properties  of  the  pulse  function  II  of  (25)  and  of  the 

triangle  function  Am  of  (26),  and  r“  of  (27)  can  be  simplified 
[5,6]  to 


Az/2 


Z  -  1  -  r| 
ran  J  47rkAz 


-Az/2 


K(z 


m+1 


Zn  ~  « 


2 [1  -  4(kAz)2]  K(z  -  z  -£) 
m  n  w 


+  ’  zn  *  »} 


(29a) 


15 


and 


[l-  IclMz]  %  P 


G  (p,  C  +  z  ;  £  +  P_)d£  d£  . 
m  n 


-Ap/2 


(29b) 


Integration  by  parts  twice  and  the  change  of  variable  z'  = z^ + £  are 

employed  to  transform  (27b)  into  (29a)  while  the  variable  changes 

z  =  z  +C  and  p'=p  +5  are  employed  to  transform  (27c)  into  (29b). 
m  n 

Also  involved  in  the  conversion  from  (27b)  to  (29a)  is  a  very  accurate 

approximation  [5,6].  By  making  use  of  properties  of  II  and  the 

variable  changes  mentioned  above,  the  expressions  of  (28)  for  Fas  , 

ran 

c 

Y  ,  and  Y  can  be  converted  to 
ran  mn 


az/  i. 

C  '  -  &  |j5  \  cs(P.o;v?)de] 


-Az/2 


Ap/2 

Y  -  j  £  [  Ga(p  ,0;p  +£)d£  , 

mn  n  j  m  n 

-Ap/2 


(30a) 


(30b) 


GC(p  ,0;p +£)d£ 
pn+t,  m  n 


-Ap/2 


(30c) 


To  facilitate  subsequent  discussion  it  is  convenient  to  express 
(27a)  and  (28a)  as  a  single  matrix  equation  of  the  form 


Zmn  rmn  J*n  0 

ras  Y  -  yc1  [v  h8C 

mn  mn  mn]  in  m 


*,»  v  ’-’V  v  v  \**\*"*v  •  •’  *  •**  »*•  .*•  .** . 

v  C  w  ■w.  C  I.'  cT  ft.*  ^  •.  r n  ft_*  ft  i  ft  . 


ts  *  4v  , 


Analytical/numerical  procedures  for  computing  the  matrix  elements  from 
(29)  and  (30)  are  presented  below. 


Analytical  Aids  in  the  Numerical  Evaluation  of  Matrix  Elements 


In  this  section  are  presented  analytical  procedures  which  greatly 

enhance  the  accuracy  and  efficiency  of  methods  for  numerical  evaluation 

of  the  matrix  elements  of  (31).  In  fact,  without  the  analysis  outlined 

below  it  is  unlikely  that  sufficient  accuracy  of  matrix  element  values 

could  be  attained  to  allow  acceptable  numerical  solutions.  Z  is 

mn 

omitted  from  the  discussion  below  since  procedures  for  its  evaluation 
can  be  found  elsewhere  [  4,6  ] . 

S3 

As  a  first  step  in  casting  T  into  a  form  suitable  for  computation 

mn 

we  point  out  that  by  making  use  of  the  analysis  presented  in  [  7 ]  one 
S3 

can  show  that  F  of  (29b)  can  be  reduced  to 
mn 


r 


sa 

mn 


Az 

» 

(1  + 

-Az 


)tgm,n+Js(C) 


gm  _  M)dC, 
m.n-’j 


(32a) 


where 


8»,nil5<C) 


-jk  R  , 
o  m.nt^ 
e _ 

R 

m.ni^s 


d<p’ 


(32b) 


and 


R2  .  *  (c  +  z  )  2  +H2  +  (p  ±  h  Ap)  2  -  2K  (p  ±  k  Ap)cos4>'  .  (32c) 

m.ni^  m  n  n 
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Evaluation  of  (32)  is  essentially  the  same  chore  as  evaluation  of 

matrix  elements  in  cylindrical  antenna/scatterer  theory  and  it  can  be 

achieved  by  methods  developed  for  that  purpose.  In  fact  for  n  =  l, 
g  ,s  g  ,  is  a  special  case  of  the  exact  kernel  for  the  circular 

cylinder  which  for  z^=  0  is  logarithmically  singular  at  £  =  0.  Hence 

from  the  simple  expression  of  (32) ,  one  can  avail  himself  of  any  of 

the  numerous  procedures  [4,6]  developed  for  cylindrical  antenna/ 

S3. 

scatterer  theory  to  evaluate  T 

mn 

3S 

A  form  of  r  suitable  for  numerical  computation  is  developed 
by  noting  that  an  interchange  of  differentiation  and  integration 
in  (30a)  is  valid  for  pfA..  Since  the  bracketed  expression  in  (30a) 

3S 

is  evaluated  only  at  points  p-p  >^+^sAp,  T  can  be  rewritten 

m —  mn 

as 


h  TT 


„as 
'  mn 


1 

8tt2 


-jkD 


m 


Kn (z '  )  (1  +  jkD^)  (pm  -  Acostf) 1 ) 


-h  -ir 


"d3" 

m 


d^’dz’ 


(33) 


h  TT 


_  1 
8tt2 


n  (z’)  i 

n  as 


-h  -TT 


d4>'dz’ 


where 

D  *  [p2+-'l2-2p  Acos  <J>'  +  z'2]*5. 
mm  m 


The  integrand  in  (33)  is  bounded  and  one  can  evaluate  the  integral 
numerically.  Because  of  the  small  size  of  the  coax  aperture, 
however,  the  observation  point  may  be  quite  near  the  cylinder  surface 
relative  to  the  length  of  a  source  segment  on  the  cylinder.  Thus 
the  integrand  of  (33)  may  be  very  highly  peaked  for  the  cylinder 
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source  segment  in  which  |z'|  <  Az/2,  To  facilitate  evaluation  of 
(33)  in  this  case  the  integral  is  expressed  in  a  form  more 
appropriate  for  numerical  calculations.  We  note  that 


Tas  kD  -H3 
m 


—  +  ~— 
D3  2D  , 
m  m' 


(p  -  H  cos  (j> ' )  =  I ' 
m  as 


(34) 


3S 

Thus  when  n=  (Ng+l)/2,  r  is  expressed  as 


„as 

mn 


Az 

2  7T 

1  f 

8  IT2  J 

-Az  — 7T 


(Ias' 


I'  )dtj>'dz' 
as 


1 

8tt2 


Az 

2 


I’ 

as 


d(J)  'dz ' 


-Az  -TC 


(35) 


The  highly  peaked  behavior  of  the  integrand  in  (33)  is  then  localized 

to  the  integrand  of  the  last  integral  in  (35).  Furthermore,  integration 

with  respect  to  the  variable  z’  can  be  performed  analytically  on  I' 

as 

to  yield 


1 

8tt2 


I'  dz’dcj)' 
as 


+  8TT2" 


f  -r2' 


as 


-IT 


d<j> ' 


(36) 


where 


,1’ 


Az[  (p  -^i)  +  4(l-cos0 1 )  ] 
m 


as 


[(p  -/i)2 +2p  4(l-cos(j)')][(p  -/t)2  +  2p  *(l-cos<t>’)  +  (~-)z]h 
mm  mm  i 


(37a) 


2  * 

Iag  =»  k2  (p^- -'tcostj)')  •{  Jin 


+  [p2  +  K1  -  2p  A-cos^'  +  (^r-)2]*5 
Z  m  m  z 


-  4  £n(p2  +  H2  -  2p  Acos<pr) 


)  . 


(37b) 


The  integrand  I  is  slowly  varying  with  respect  to  <J> '  and  can  be 

3.S 

X  ? 

evaluated  numerically.  The  integrand  I  ,  however,  is  still  highly 

as 

peaked  near  We  therefore  subtract  the  small  argument  behavior 

1 ' 

of  I  under  the  integral  sign  and  add  it  to  the  expression  for 
as 

3S  2l  S 

r  as  a  separate  integral,  so  that  F  is  finally  given  by 
mn  mn 


Az  IT 


8^  j  2  I  (Ias-IL>d*’dz’ 


-Az  -7T 
2 


if  1’  1"  if?' 

+  8^  W*' 


+  8^2  f  Iasd41’  ."-Ws+1)/2 


1 1  2 ' 

where  I  ,  I'  ,  I  ,  and  I  are  defined  in  (33) -(37)  and  where 
as  as  as  as 


nil  Az[(p  -fl)  +  i  X<p'Z] 

_  _ m _ Z _ 

as  t(p  -n)2  +  &)2}\(q  -4) 2  +  p  *f>'2] 
m  z  m  m 


(39a) 


i  h  a  p  a 

I  d(j)f  =  - - - T - 7 —  7T  +  —  ■  — 

p.[(p  -V2  +  Joji 


P  +  K  TT/p  Jt 

tt  +  —2 — _  tan'1  - SL_ 

/pT  Pm"  ^ 

m  v 


(39b) 


For  nj*  (N  +l)/2,  of  course,  T  may  be  evaluated  directly  from  (33) 
s  mn 


via  numerical  integration. 


sty 


In  view  of  the  above,  Y  of  (30b)  becomes 

mn 


Ap/2  f  TT 


Y  =  j  — 
mn  ryrr 


-Ap/2  0 


+Tv^K(6n„>h- 

v  mn  m  n  i 


When  m  =  n,  K  is  unbounded  for  t,  0  but,  fortunately,  the  small  argument 
form  of  K  can  be  integrated.  That  is. 


(2pn+C)  nn 


K(8  )  -r—t  -  —  In 

nn  £  -+  0  p 


t  ’-n8rd5 

n  -Ap/2  n-» 


Hence,  Y  is 
nn 


Y  *  JL  M  i  _  £n  [_^£ 

nn  JnTr  I  P  16p 

n  n 


-ao/2 


(2THTK<6nn>+f  « 

n  n  n 


Ap/2  7T r 
f  f  cos^1  e 

J  i  R 

-Ap/2  0 


dc^'d^ 


The  unbounded  part  of  the  elliptic  integral  is  added  to  and  subtracted 
from  the  integrand  of  (42) .  The  integral  of  the  added  term  becomes 
the  first  term  of  (45)  while  the  subtracted  term  remains  in  the  integrand 
of  the  first  integral  of  (45).  Therefore,  both  integrals  of  (45)  have 


very  slowly  varying  integrands  and  can  be  integrated  numerically  with  ease. 


The  integration  indicated  in  (30c)  can  be  performed  analytically 


and,  as  a  result. 


YC  becomes 
mn 


mn  .bp 

n  Jin  —  m 
c  a 


~~~  In 


n+*j 


n-% 


1  -  T  e 


-j2kcL 


1  +  r  e 


-j2k  L 
c 


-  —  l  — L  4- $  cp  )f*  (p 

nc  q=l  XqWq  dp  ^  1 


wvw]  (46) 


in  which  p  =  a  +  nAp  and  p  ,  =  a  +  Ap(n-l)  and  in  which  -r-  $  (p  ) 
n+*5  n-*s  dp  q  m 


means 


Ip  V’> 


P=P 


m 


The  series  in  Y  must  be  summed  efficiently  if  the  numerical 
mn 

procedure  described  here  is  to  be  practicable.  This  is  particularly 

true  in  regard  to  Y  because  of  the  need  to  solve  a  transcendental 

mn 

equation  to  determine  the  transverse  eigenvalue  kfc^  of  (6)  prior 
to  the  computation  of  each  series  term.  To  this  end  Rummer's 
transformation  [9]  is  used  to  accelerate  convergence.  To  employ 
this  technique,  one  determines  the  q*"*1  term  for  large  q  of  his 
series  and  then  finds  another  series  with  known  sum  and  the  same 
large- index  q  term.  Even  though  Rummer's  transformation  is  applied 

to  the  series  of  tne  present  problem  under  the  condition  that  the 
guide  is  below  cutoff  of  the  higher-order  modes,  k2<k2^,  which  is 
the  case  of  major  practical  interest,  the  procedure  is  immediately 
extendable  to  the  cases  in  which  higher-order  modes  do  propagate. 


One  can  show  readily  that 


and  subsequently  that 


X  N2  dp  %^Pm^  ^VPn±^  q-«°  Cmn±*s 


vpm  pn±*s  /(a+Ap[m-*5] )  (a+Ap[  (n-Jj)!^] ) 


9mn±»s  "  W  (pm  +  Pn±is  “  2a)  "  W  C  CnH-n-1 ) ±*5]  , 


mn±*5 


In  [10]  the  series  below  and  its  closed-form  sum  are  listed: 


where 


£  S  (a)  =  4(a),  0<a<2ir 


„  ,  N  sin  qa 
Sq(a> 


*  J^-  -  ?,n  (2  sin  -^)j  s inn ,  0<ci<2tt 


Notice  that,  apart  from  the  factor  C^,  exhibits  the  form  of  the 
right  side  of  (47b)  when  q-+°°.  Applying  Rummer' s  transformation  to  the 


series  in  Y  ,  one  arrives  at 


<■  Mn-V  \l  +  r  e  J  C_J 


H  In—  m 
c  a 


which  converges  in  very  few  terms.  Absolute  values  are  used  above 
because  t ^  can  be  negative  and,  as  indicated,  (49c)  is  valid  only 

for  0<a<27T.  Since,  for  a<0,  sin  a  ■  -sin  [a|,  we  can  write 


sin  a  ■  -2—  sin  la  I 


for  all  a 


( 


This  can  be  used  to  advantage  to  rewrite  (49)  as 


y  sin  q  a  a  “  sin  q  laI  a  1  „  , 

1  - r  *  -  >  - 5— ; -  =  -  7 —  in  (2  sin 

q-2  q  -1  |a|  q-2  q  _1  lal  4 


for  -2Tr<a<27T 


l  S  (a)  =  — ~  -  £n(2  sin  sin  lal,  -2ir<a<2ir 

q  I  -  I  4  2  '  ' 
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4.  SAMPLE  NUMERICAL  RESULTS 

In  this  section  are  included  sample  numerical  results  obtained 
from  a  computer  program  implementing  the  techniques  developed  in 
Section  2.  In  Figs.  4,  5,  and  6  are  shown  current  distributions  and 
aperture  field  distributions  for  coax-fed  monopoles  above  ground  plane 
with  the  shape  factors  ft  =  8.54,  9.34,  and  9.92,  respectively,  where 
ft  *  2  2n  (2h/A.) .  The  monopoles  reside  in  air  and  have  radii  fi  = 

3.175  mm,  and  equal  to  the  radius  a  of  the  inner  conductor  of  the  air- 
filled  coaxial  line  with  b/a  =  3.  The  operating  frequency  is 
f  =  663  MHz.  The  computed  current  distributions,  which  are 
normalized  in  each  case  to  the  voltage  in  the  aperture  V  ,  are 
compared  with  the  experimental  data  (shown  by  dots)  obtained  by 
Mack  [11].  The  agreement  is  seen  to  be  very  good.  Although  no 
measured  data  are  available  to  corroborate  the  computed  aperture 
field  distributions,  it  is  seen  that  they  exhibit  the  expected 
behavior  near  the  inner  and  outer  conductors  of  the  coax. 

In  Figs.  7  and  8  are  shown  computed  results  for  a  very  thick 
monopole  (ft  =  2.77)  immersed  in  fresh  water  (e^  =81,  o  -  0)  and 
sea  water  (e^  =  81,  a  =  4  S/m),  respectively,  and  driven  by  an  air- 
filled  coaxial  line.  The  other  parameters  are  A.  *  a  =  2.78  cm,  b/a  = 
and  f  *  300  MHz.  It  is  observed  that  in  the  lossy  water  case  the 
monopole  current  decays  rapidly  away  from  the  driving  point.  Also, 
the  driving-point  admittance  is  much  higher  in  this  case  than  in  the 
case  of  fresh  water. 
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vCvivlv 


\  ■*. 


I  /Va  [mA/v] 

(a) 


Current  and  aperture  electric  field  for  one-quarter  wavelength 

coax-fed  monopole  in  air  ( - computed,  ...measured),  (a)  mono 

pole  current,  (b)  aperture  electric  field. 


I  /V  [mA/V] 


(a) 

Current  and  aperture  electric  field  for  three-eights  wavelength 
coax-fed  monopole  in  air  ( - computed, .. .measured) .  (a)  mono¬ 

pole  current,  (b)  aperture  electric  field. 


I  /V  [mA/V] 
(a) 


Current  and  aperture  electric  field  for  one-half  wavelength 

coax-fed  monopole  in  air  ( - computed,  ...measured).  (a)  mono 

pole  current,  (b)  aperture  electric  field. 


I/V  [tnA/V] 
(a) 


Current  and  aperture  electric  field  for  coax-fed  monopole  in 
tap  water.  (a)  monopole  current,  (b)  aperture  electric  field 


I/V  [mA/V] 
(a) 


Current  and  aperture  electric  field  for  coax-fed  monopole  in 
lossy  medium.  (a)  monopole  current,  (b)  aperture  electric  field 


Finally,  in  Fig.  9  is  shown  a  plot  versus  frequency  of  the 


driving-point  admittance  (normalized  to  Yq  =  0.029  S)  of  a  coax-fed 
monopole  of  radius  1  =  a  =  0.31  cm  and  height  h  =  1.65  cm  immersed 
in  tap  water  (e  =  82,  o  =  0.11  S/m).*  The  driving  coaxial  line  has 
b/a  =  2.29  and  is  air-filled  except  for  a  teflon  bead  (e  *  2.1) 
of  length  L  =  0.99  cm  at  the  ground  plane.  The  computed  results 
(shown  by  boxes)  are  compared  with  experimental  data  (shown  by 
pluses)  obtained  by  Harrison  and  Butler  [12].  Although  the  computed 
conductance  is  generally  larger  than  that  determined  experimentally, 
there  is  excellent  agreement  at  peak  values.  A  major  difference  of 
some  concern  is  that  the  computed  frequency  for  zero  phase  angle 
does  not  coincide  with  the  computed  frequency  for  peak  conductance. 

In  the  experiment,  peak  conductance  and  zero  phase  angle  occurred 
at  approximately  the  same  frequency.  The  wavelength  in  the  water 
at  500  MHz  is  approximately  2.6  inches,  and  the  monopole  height  is 
approximately  0.65  inch  while  its  radius  is  0.25  inch.  Fabrication 
of  an  antenna  of  these  dimensions  to  close  tolerances  relative  to 
the  above-mentioned  wavelength  is  difficult,  especially  so  since 
the  monopole  must  be  suitable  for  underwater  operation.  Construction 
inaccuracy  is  likely  a  contributor  to  the  failure  of  computed  and 
measured  data  to  agree  better. 


Actually  and  0  are  slightly  frequency  dependent  over  the  frequency 

range  under  consideration;  these  representative  values  are  used  in 
the  computations  for  convenience. 


Frequency  (GHz) 


Normalized  (with  respect  to  YQ(-0.029  S))  admittance 
of  coax-fed  monopole  in  water. 
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